Exact relativistic models of thin disks around static black holes in a magnetic field 
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Using the Rayleigh criterion we derivate for circular orbits the stability conditions of the particles 
of the streams. 
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I. INTRODUCTION 

There is a strong observational evidence that active galactic nuclei ( AGN) , X-ray transients and gamma-ray bursts 
(GRBs) are associated with black holes that accrete matter via a surrounding accretion disk. The exact mechanism 
by which these phenomena are produced involves the interaction between a rotating black hole, the accretion disk 
and the electromagnetic field [ll-Q- However, magnetic fields play a key role in understanding of these process. The 
enormous observed energy in AGN is related with the presence of magnetic field in these nuclei ■ The existence 
of radio jets is also attributed to the presence of strong magnetic fields in centers of AGN and quasars 

A general exact relativistic model that describes such astrophysical objects require an exact solution of coupled 
Einstein-Maxwell field equations that represent the superposition of a Kerr black hole with a stationary disk and 
electromagnetic fields. As a first approximation one could consider a static system composite by a Schwarzschild 
black hole and a thin disk immersed in a magnetic field. Exact solutions of the Einstein equations representing the 
field of a static thin disks without radial pressure were first studied by Bonnor and Sackfield ,13], and Morgan and 



Morgan [14{ , and with radial pressure by Morgan and Morgan [15| . Several classes of exact solutions of the Einstein 
field equations corresponding to static thin disks with or without radial pressure have been obtained by different 
authors TB;22|. Rotating thin disks that can be considered as a source of a Kerr metric were presented by Bicak and 



Ledvinka [23|. while rotating disks with heat flow were were studied by Gonzalez and Letelier [24j. Also thin disks 
in presence of electromagnetic field have been discussed as sources for Kerr-Newman fields [25j . conformastationary 
metrics [26[ and for magnetostatic axisymmetric fields in [27i - [29| . The static superposition of a disk and a black hole 
was first considered by Lemos and Letelier (30l - [33 |. 

In this work we consider the exact static superposition of a Schwarzschild black hole and a thin disk in presence 
of a magnetic field. The method used to include the magnetic field is the well-known complex potential formalism 
proposed by Ernst [HI, [34| , using as seed solutions simple vacuum spacetimes representing the field of a thin disk and 
a black hole. 

The paper is organized as follows. In Sec. II we present a summary of the procedure to obtain models of thin 
disks with a purely azimuthal pressure and currents for the Einstein-Maxwell equations. In order to have a stable 
configuration in absence of radial pressure, the matter in the disks also in interpreted as made of two pressureless 
(dust) streams of counterrotating charged particles (counterrotating model) moving along electrogeodesic. Using the 
Rayleigh criterion we derivate for circular orbits the stability conditions of the particles of both streams. 

In Sec. Ill we consider disk constructed from magnetized Weyl solutions. We shows that if the gravitational field 
is continuous across the plane of disk but its normal derivative is discontinuous, the magnetized solutions constructed 
using the Ernst' method can be interpreted as the exterior gravitational field produced by a dislike source immersed 
in a magnetic field. 

In Sec. IV we present the magnetized version of the Schwarzchild solution in prolate coordinates. In Sec. V the 
formalism for superposing disks and black hole in the vacuum [3lj is extend to the case of magnetized Weyl solutions. 
In Sees. VI- VII we analysis two models of disk solutions with black hole in presence of magnetic field. The first model 
is based in the Kuzmin-Chazy-Curzon inifnite disk and the other in the first Morgan-Morgan finite disk. Finally, in 
Sec. VIII we summarize and discuss the results obtained. 

II. ELECTRO- VACUUM RELATIVISTIC THIN DISK 

The simplest metric to describe the exterior gravitational field produced by a static axially symmetric body is the 
Weyl's line element [43[ 

ds 2 = - e 2 ^dt 2 + e- 2 ^[p 2 dip 2 + e 2A (dp 2 + dz 2 )}, (1) 

where (t, tp, p, z) are the Weyl canonical coordinates, and i\> and A are functions of the coordinates p and z only. For 
the coordinates we also use the notation (x , x , x 2 , x 3 ) = (t, ip, p, z). The vacuum Einstein-Maxwell equations, in 
geometrized units such that G = c = 1, are given by 

Rab = S7TT ab , (2a) 

(2b) 

F ab , b = 0, (2c) 
F ab = A b>a - A a<b , (2d) 



Tab — ~r- 
4w 



F ac F b c — -^g a bF c dF c 
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where T a b is the electromagnetic energy-momentum tensor, F a b is the electromagnetic field tensor, and A a = ((f), A, 0, 0) 
is the four-potential, where <fi is the electric potential and A the magnetic potential which are also functions of r and 
z only. Further, ( ). a = d/dx a , and ( ) ;a means covariant derivate. 

For the metric (JXJ) and in magnetostatic case, the Einstein-Maxwell equations are 

V • [p- 2 fVA] = 0, (3a) 
/V 2 /-V/-V/ + 2p- 2 / 3 VA-VA, (3b) 

K = P (V> 2 ,p ) + ~ (A 2 1P -A\ z ) /, (3c) 
2 

A, z = 2p^, P i>, z +-A, p A, z f, (3d) 

where / = e 2 ^ . The equations (|3a|) y (|3bl) are equivalent to [H, [34| 

fA£ = (V£ + 2$*V$) ■ V£, (4a) 

/A$ = (V£ + 2$*V$) ■ V$, (4b) 

where A and V are the standard differential operators in cylindrical coordinates, / = e 2 ^, and £ and $ are complex 
potentials which in the case static £ = £ * and for the case magnetostatic $* = — <f> . The above equations are called 
Ernst equations. The metric functions are obtained via 

/ = £ + $$*, (5a) 



where V2~C — P + so that v2t?£ = 9. p — z<9. z , and the magnetic potential A is related to $ via 

A, c = ij(hn$) x . (6) 

In order to obtain a solution of the Einstein-Maxwell equations (|2"aj) - (f2"c|) representing a thin disk at z = 
with electric current, we assume that the components of the metric tensor and of the electromagnetic potential are 
continuous across the disk, but have discontinuous first derivatives in the direction normal to the disk, which can be 
cast as 

bab = [g a b,z] = 9ab,z\ z=Q+ - 9ab,z\ z=0 _ = 2 9ab,z\ z=0+ , (?) 

a b =[A b , z ]=Ab, z \ z=0+ - A h , z \^ _ = 2A M |„ 0+ . (8) 



By using the distributional approach [44M46| or the junction conditions on the extrinsic curvature of thin shells [47, 48], 
the Einstein-Maxwell equations yield an energy-momentum tensor and an electric current density 

Tab = T Q r + = T a r + Q ab S(z), (9a) 
Ja = jaS(z) = -g zz a a S(z) = -2e 2 ^>-^A a , z 8{z), (9b) 

where S(z) is the usual Dirac function with support on the disk, T^ b m is the electromagnetic tensor defined in Eq. 
(I2b|) . and 

Qb = lib az S z b - b zz 5 a b + g az bi - g* z b a b + b c c (g zz S a b - g az 5 z b )} (10) 

is the distributional energy-momentum tensor. j a is the electric current density on the disk. The "true" surface energy- 
momentum tensor (SEMT) of the disk, S ab , and the "true" surface current density, j a , can be obtained through the 
relations 

Sab = J TX 1 d Sn = e A ~* Qab , (Ha) 
j a = I Ja ds n = e A ~^ Ja , (lib) 
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where ds n = y/g Z z dz is the "physical measure" of length in the direction normal to the disk. For the metric ([T]), the 
nonzero components of are 

Sg= 2e^- A {A, z - 2^}, (12a) 

Sl= 2e^ A A, z , (12b) 
and the only nonzero component of the current density j a in the magnetostatic case is 

h = -^~ A Az, (13) 

Z7T 

where all the quantities are evaluated at z = + . 

In order to give physical significance to the components of the energy-momentum tensor S b and the electric current 
density j a we project them onto the orthonormal tetrad e^ a ) b = {V b , W b , X b , Y b }, where 

V a =e^ (1,0,0,0), (14a) 

W a = — (0,1,0,0), (14b) 
P 

X a = e^ A (0, 0,1,0), (14c) 

Y a = e* _A (0, 0,0,1). (14d) 

In terms of this tetrad (or observer with four-velocity V a ) the surface energy density e, the azimuthal pressure p v of 
the disk are given by 

e = S\, V v = S$, (15) 

and the azimuthal current density j 

j=W%. (16) 

Thus we have a disk only with pressure and electric current in azimuthal direction. Because there is no radial 
pressure or tension to support the gravitational attraction, the matter distribution is unstable. In addition, since the 
spacetime is static we have no rotation. In order to have a stable configuration in absence of radial pressure, we need 
assume the counterrotating hypothesis, that is the matter in the disk is considered made of two pressureless streams 
of counterrotating charged particles, i.e., that circulate in opposite directions. Even though this interpretation can be 
seen as merely theoretical, there are observational evidence of counterrotating matter components in certain types of 
galaxies [35M39| . We assume 

gab = gab + gab ^ ( 1?a) 
f=f++f-, (17b) 

where 

S^ b = e± t4t4, (18a) 
f ± =cr±u a ± . (18b) 

e± are the matter densities of each stream, a± the electric charge densities, and u± the normalized four-velocities 
(u±u a ± = —1), which for circular orbit are (u±) = (u±,u±,0, 0) = u±(l,w±, 0, 0), where uj± = v}±/u\. are the angular 
velocities of each stream. Furthermore, we can assume in the disk the continuity equation (the Bianchi identity) 

Sf. b = F\J b ± . (19) 

From Eqs. (|18a[) . (|18b|) we have e±u±u±. b = <J±F a b u b ± ) i.e., each stream follow a electrogeodesic motion. For circular 
orbits, the equation for the electrogeodesic motion of the particles in the disk is given by 

^e±3a6,pw±w± = -<J±F pa u a ±1 (20) 
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and in the case magnetostatic reads 



~e±u°±{g VVtP u 2 ±+ g tt ,p) 



-a±A tP u)±, 



where 



u± obtains normalizing u±, that is 



(4) 2 = - 



i 



g w u± + g t t 



Using (ITTal) - (|18bj) and (gT]) we obtain 



M^ 1 



e 



1 



With respect to the orthonormal tetrad (TT4")) the 3-velocity has components 



(») 



and for equatorial circular orbits the only nonvanishing velocity components is given by 



SW. ,2 

w ± 

gu 



(21) 
(22) 

(23a) 
(23b) 
(23c) 

(24) 
(25) 



and represents the circular speed of the particles as seen by an observer at infinity. In fact, when p — ¥ oo, ip — > 
and t> 2 = p 2 uj 2 = p-jS where </> is the Newtonian gravitational potential, and v is the circular speed, which represents 
the speed of a test particle in a circular orbit at radius p. Using (|23a[) one finds that the speed v of counterrotation 
(rotation curves or rotation profile) of the particles in the disk is given by 



e 



(26) 



To analyze the stability of the particles of the two streams in the case of circular orbits in the equatorial plane we 
use an extension of Rayleigh criteria of stability of a fluid at rest in a gravitational field [40 42] . The method works as 
follows. Any small element of the matter distribution analyzed (in our case a test particle of the streams) is displaced 
slightly from its path. As a result of this displacement, forces appear which act on the displaced matter element. If 
the matter distribution is stable, these forces must tend to return the element to its original position. 

The relativistic Lagrangian for a test particle in presence of a gravitational and magnetic field is given by 



£± = -g a bu a ±u b ± + a±Au^, 



where a± is the specific electric charge. For magnetized Weyl fields we have two constants of motion 

E± = -gttu±, 

L± = g vv w±u\ + a±A, 



(27) 



(28) 
(29) 



where E± represents the relativistic specific energy and L± the specific angular momentum. The motion equation 
(|2T|) can be cast as a balance equation 



f p gtt, P E 2 ± 



eg pp A, p lo±u° ± 



°g vv , P (L± - cr±A) 2 
2flL, 



(30) 



where the term first on the left-hand side represents the gravitational force F g , the term second the the Lorentz force 
Fl, and the term on the right-hand side the centrifugal force F c (p) = F{p 1 L±{p)) acting on the test particle. So 
we have a balance between the total force F(p) = F g + and the centrifugal force. We now consider the particle 
to be initially in a circular orbit with radius p — p$ and we slightly displace it to a higher orbit p > pq. The 
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angular momentum of particle remains equal to its initial value L±q — L±(po) which implies that the centrifugal 
force in its new position is F c (p, L± ). In order that the particle returns to it initial position must be met that 
F(p) > F c (p,L± ), but according to the balance equation (f30|) F(p) — F c (p,L±) so that F c (p,L±) > F c (p,L± ), and 
hence (L± — <J±A) 2 > (L±o — &±A) 2 . Using the expression for L± (|29p and defining la function h± — gip V uJ±u±, 
follows that h±(p) 2 > h±(po) 2 - The quantity h± can be written as 

*± = 7?=^ (31) 

and has the same form that the specific angular momentum in the vacuum, being the true specific angular momentum 
corresponding to expression ([29)) . By doing a Taylor expansion of h±(p) around p = po one finds that the condition 
of stability for equatorial circular orbits is 

h±h ± , p > 0, (32) 

or, in other words, h 2 ± > 0. Thus when the counterrotating hypothesis is assumed, the stability of the disks is 
equivalent to the stability of the particles of the two streams. 



III. DISKS FROM MAGNETIZED WEYL SOLUTIONS 



A class of exact solutions to the Einstein-Maxwell equations can be obtained taking £ as function of <£>. From l]4ap 
and (I4b[) it follows that 

(R/£ + |$| 2 )0V<&-V4 = O, (33) 

and hence that £ is a lineal function of $. Using the boundary conditions £ — > 1 and $ — > at infinity, we obtain 

£ = l-2q- 1 $ 1 (34) 
where q is a complex constant. With the change of variable 

f = fri' (35) 



then (j34|) implies 



*=«TT)' (36) 

and the Ernst equations read 

Kr-(i-<7<7*)]v 2 e = 2rve-ve (37) 

Then making £ = (1 — qq*) 1 ^ 2 ^ the equation (|3"T)) takes of form 

(||* - 1)V 2 | = 2|*V| • V|, (38) 

which is the Ernst equation in the vacuum [34| . So given a solution of the Einstein field equations in vacuum | (seed 
solution) we can construct a solution of the Einstein-Maxwell field equations. A solutions for this equations is 

|=-e ia cothV;, (39) 

where the function tp satisfies the Laplace's equation 

V 2 ^ = 0. (40) 

The case a — corresponds to the well-known Weyl vacuum solutions and thus the magnetized solutions (taking 
$ imaginary) built from them can be called magnetized Weyl solutions. The metric functions and magnetic potential 



7 



are given by [2£ 



/ = 7 7f , ( 41a ) 

{fi + 1) er§ - (/3 - 1) 

A = A 



{(^ 2 , p -V> 2 , z )dp + 2^, p ^dz}, (41b) 
p{—i/>, z dp + V>, p dz} , (41c) 



where /? = s/l + 6 2 , being 6 is the parameter that controls the magnetic field, and A is the metric potential A 
corresponding to the seed solution, that is taking £ = e 2 ^ and $ = 0. 

We consider magnetized Weyl solutions representing the field of a disk. The components of the momentum-energy 
tensor of the disk Q h a can be obtained by integration of the field equations (|2al) writing T a t, as (l§a|) [l7j 

z— 0+ f-z— 0+ />z=0+ 

i? a6 dz = 8tt( / T a e ' m rfz + / Q ab S(z)dz). (42) 



The nonzero components of T°£, m are 



z=0- Jz=0- Jz=0- 

clm 



e 2(3V-A) 

8np 2 



2S3r = -5z^-(^+^). ( 43a ) 



Tj}" 1 = -L e 2 W- A )(A 2 + A 2 J, (43b) 

87T ,P 



^ m = -^ m =^(Af p -^), (43c) 

T2r = i^ A >o A >- (43d) 

T elm = ^ e 2(3^A) ( ^ + ^ )) (44a) 

T i r = ^pV^- A )(^ + ^) ! (44b) 
r 2 1 m = -T 3 1 m = ^e 2 ^ - ^ 2 p ) (44c) 

T 23 m = -|^e 2 ^>,*. (44d) 

Now the continuity of the metric functions across the disk implies that they are even functions of z and the 
discontinuity of its first derivatives in the direction normal to the disk means that these are odd functions of z. 
Thus ip >z is odd and its square is even and hence continuous across the disk. It also follows that ip,p is continuous. 
Accordingly the terms 

z— 0+ f-z— 0+ f-z— 0+ f-z— 0+ 

T^dz = / T^dz = / T^dz = / T^dz = (45) 

z=0- Jz=0- Jz=0- J z=0- 

are zero due to the continuity of the metric, ■0 2 p and ip 2 z across the disk. For the other component we have 

T^dz = ~—e 2 U, P / t*dz = - — e 2 %M = 0, (46) 

which also vanishes because continuity of the metric and ip,p across the disk. Therefore, the first term on right hand 
side of (|42|) is zero and in consequence Q b a has the same form that the vacuum (fT0|) . 

On the other hand, on the disk the Maxwell equations dbF ab = AirJ a , where 'bar' denotes multiplication by y/— g, 
are given by 

- An] v S(z) = 8,(0" A,,) + d p (g">A, p ), (47) 



2e 



and using (|41cl) we have 
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but again using (|41cl) we have 



Integrating through the disk we have 



4xj v 6{z) = bpd z (g zz ip, P ) - bd p (pg^,z)- (48) 



- 4ttj v = bp f + d z (g zz ^, P )dz -b f + d p {pg pp ^, z )dz 

= bp(g zz i>,p)\Zt - bd P (pgP p f ° + hdz) 
= bp[g zz i P )}-bd P (pg pp [^}) 

= 0, (49) 

where the terms on the right hand side vanish due to the continuity of the metric and ij> iP across the disk. Since 
electric current density is zero on the disk the source of the magnetic field is non planar but of a different origin 
such as a remnants or fossil magnetic field (52j . or can come from external sources, such as the presence of a nearby 
magnetars or neutron stars. Thus, if we have a continuous gravitational field across plane z — but its normal 
derivative is discontinuous, the magnetized solutions constructed using the Ernst' method can be interpreted as the 
exterior gravitational field produced by a dislike source immersed in a magnetic field. 

For completeness, let us also analysis the electrostatic case, i.e., for $ real. The electric potential <f) is given by 



(1 + o)e-^ + (1 - a)ev 



(50) 



(51) 



me where a = y/l + p 2 , being p the parameter that controls the electric field, but as 

, = 4p^ 

[(l+a)e-^ + (l-a)e^] 2 ' 

then the electric potential 4> is discontinuous and the current density in the disk is nonzero 

.it = -^~ A ^ (52) 
whereas in terms of the tetrad (|140 the electric charge density a is given by 

a = -V°j . (53) 
Therefore, in this case the solutions can be interpreted as the gravitational field produced by a charged disk. 



IV. MAGNETIZED SCHARWZCHILD SOLUTION 



In terms of the prolate spheroidal coordinates (x, y) which are related to Weyl coordinates (p, z) by 

p 2 = k 2 (x 2 - 1)(1 - y 2 ), z^kxy, k = const (54a) 



2kx = r+ + r_, Iky = ?*+ — r_ 



p 2 + (z±k) 2 , 



(54b) 



with x > 1 and — 1 < y < 1, the asymptotically flat general solution of Laplace's equation (|40| can be written as 

oo 

j> = C 2n P2n(y)Q2n{x) , (55) 



n=0 



where C2„ are constants, P2n(y) are the usual Legendre polynomials, and Q2n{ x ) the Legendre functions of the second 
kind. The function metric A can be obtained by integrating 



A.* = 



\v = 



i-y 2 

x 2 — y 2 
x 2 - 1 
x 2 — y 2 



x(x 2 - - x(l ~ y 2 )f y - 2y(x 2 ~ l)txtv 
y(x 2 - l)^ x - y(l - y 2 )4> 2 y + 2x(l - y 2 )tJ, v 



(56a) 
(56b) 
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For magnetized solutions constructed from the Ernst's method the magnetic potential is given by |4£ 

oo 

C2r, 



-kb (x 2 



l ) E 2^+i Q ' 2n ^ [yP2n[y) ~ P2 "-^' 



(57) 



with Q' 2n {x) = dQ2n{x)/dx. 
For n = and cq — 1 we have 



and therefore 



n=0 



v 2 \x + l 



,i> - 



2(x 2 - I) 1 ' 2 



{fi + l){x + l)Va - (/? - i)( x - 1)V2 



A = - In 
2 



x l - 1 



— y 2 



A = fc6(y + 1). 



(58) 



(59a) 

(59b) 
(59c) 



In the absence of magnetic field b = this spacetime is the Scharwzchild metric |50| which represents the external 
gravitational field produced by a static black hole of mass rh and parameter k = rh. Thus the case b ^ correspond 
to a magnetized Scharwzchild solution. The parameter k is related to the mass of the source as follows. For the line 
element (JTJ, in terms of the spherical coordinates, 



ds 2 = - e 2 ^dt 2 



e- 2 ^\r 2 sin 2 



+ e 2A (dr 2 +r 2 d9 2 )} 



the total mass 15l| 



in = — - lim / K a ' b dS ab 

47T ' 



(60) 



(61) 



(where S r is a coordinate sphere with radius r and K = dt is the static Killing vector) for an asymptotically flat 
solution is given by 



From (|ITa|) 



where 



F = 



m = lim (r 2 %p :r ). 
{(3 + l) e -* + [fi - l)e^ 



(62) 



(63) 



(64) 



(65) 



(/3 + l)e-^- (/3- l)e^ 
and using the fact that ip — > for r — ¥ oo we obtain 

m = lim (r 2 Fip r ) = lim (i* 1 ) lim (r 2 ijj, r ) — f3rh. 

Thus, since for a black hole k = rh, it follows that fc = m//3 in the case of the above magnetized solution. 

V. SUPERPOSITION OF A BLACK HOLE AND A THIN DISK IN A MAGNETIC FIELD 



We considerer the superposition of a Schwarzschild black hole (tpSi A-s) with a thin disk (ipD, Ad) in presence of 
a magnetic field. Since the metric function tp satisfies the Laplace's equation and it is linear, then the superposition 
"0 = $s + V'-D is a l so solution. The other metric function A = A is nonlinear but holds the relation 31] 



A = Ag + An + A-sd, 



(66) 
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where 

AsiJ = 2 J p{{'4)S,pi'D,p ~ 4>S,z4>D, z )dp + (tpS,p^D,z + ~<Ps,z^D, P )dz}. (67) 

Call As the magnetic potential associated to the magnetized Scharwzchild solution (|59c|) and Ad the magnetic 
potential for any magnetized disk solution, according to (|41c|) the magnetic potential for the composite system is 
given by 

A = A S + A D . (68) 

In prolate coordinates (x, y) the disk is located in y = 0, and from relation between oblate and prolate coordinates, 
we see that across the disk y has the same behavior that v, i.e., it changes of sign but does not change in absolute 
value, whereas the coordinate x is continuous. This implies that for a black hole both metric functions and its normal 
derivatives given by 

^ = h^wy (69) 

A,- 2 ^ 2 ~ 1 ) ; (70) 
k{x 2 — y 2 ) 2 

are continuous across the disk, i.e., [ips.z] — [As, z ] = 0. For the disk have that [if>D,z] = 2*1>d,z, and [Ad, 2 ] = 2Ad. z - 
With these considerations, we compute for the metric (TTJ) the nonzero components of b a b for this distribution of matter. 
For btt we have 

&*t = [<?«,*]= -2e 2 <^, z ]. (71) 

But from (|41aj) ip iZ = Fip yZ so that 

[tz] = [HA = F[i), z ] = F{[i)sA + {4>d,z}} = 2F$ Ajr . (72) 



So we obtain 



In the same way, 



In turn, 



ba = -4e 2 ^Fi> DtZ . (73) 
b vv = [g w A = -2p 2 e^[ip, z ] = -Ap 2 e~ 2 ^F^ z . (74) 



bpp = b zz = [g zz>z ] = 2e 2 ^\[A, z] - fofrj) = 2e 2 ( A -*>([A, z] - 2Fj> D ,.). (75) 
Using dnU) and (|41h>|) we have 

[A,J = [A s ,z] + [AdA + [Asu,»] 
= 2A DiZ + [A S dA 

= 4pi>D,pi>D, z + [A S dA, (76) 
but from (|67]) A S d,z = ^p{^s, P ^D,z + $s,Ad, p ), then 

[^sdA = ^P^sA^dA = 4p-0s, p -0D, 2 , (77) 

and hence 

[A,J=4p(&j (78) 

Thus 

bpp = b zz = 4e 2( - A ~^{2p^j s + i> D ),„ - F}^ DtZ . (79) 
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With the components of b a b (|73j) . (|74j) . and ([79)1 and using ([15)) , one finds for the combined system 

e = 4e^- A (F-p(^ + ^), p )^, ( 80 ) 
V V = ^- A p($ s + #d), p $d„ (81) 

where all the quantities are evaluated at z — + . 

Similarly, for the magnetic potential we have [A lZ ] — [Ag t g] + [Ad,z], but using (|41c|) [As, z ] = [bpips,p] = 0, so that 
for thin disk solutions in which Aj^ z is discontinuous the current density is 

U = -^- A A D , z . (82) 

VI. BLACK HOLE SURROUNDED BY KUZMIN-CHAZY-CURZON DISKS IN A MAGNETIC FIELD 

In terms of spherical coordinates (r,9) which are related to Weyl coordinates (p, z) by 

p = rsin6*, z — r cos 9, (83) 
with < r < oo and — it < 9 < it, the asymptotically fiat general solution of Laplace's equation (|40l) can be written 



C-lriPln (COS 9) 

111 ( 84 ) 



n=0 



where C2n are constants and P2n(cos9) are the Legendre polynomials. For magnetized solutions constructed from the 
Ernst's method the magnetic potential is given by 



A = b J2 [cos8P 2n (cos9) - P 2 „_ 1 (cos0)] . (85) 



71 = 

For n = and cq = M 

and in terms of the total mass (|65p we obtain 

> = 



M 

</> = , (86) 

r 



(87a) 



08 + l)e( M //3)/r _ (0 _ l) e -(Af/0)/r ' 



A=JWf.*n>6 t (87b) 
A = b(M//3)(cos6 + 1). (87c) 

In the absence of magnetic field this spacetime corresponds to the Chazy-Curzon solution [55|, [56|. Thus, for b ^ 
this metric corresponds to a magnetized Chazy-Curzon solution. 

Exact solutions which represent the field of a disk can be obtained using the well known "displace, cut and reflect" 
method that was first used by Kuzmin [53j | and Toomre [54| to constructed Newtonian models of disks, and later 
extended to general relativity [2(| [2l|, |23|. |24| . Given a solution of the Einstein-Maxwell equation, this procedure 
is mathematically equivalent to apply the transformation z — > \z\ + zq, with zo constant, on that solution, and the 
resulting disks are essentially of infinite extension. This method applied to the above magnetized solution produces 
a magnetized Kuzmin-Chazy-Curzon disk with 

M 

= ; (88) 

This disk is the relativistic generalization of the Kuzmin Newtonian disk in presence of a magnetic field. 
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Now we consider the superposition of a black hole and this disk. The interaction term Asd (|67|) between the disk 
and the black hole can be calculated in prolate coordinates from (|56b[) 



A 



SD 



X 2 - 1 

_ l \x 2 — y 2 



j {[y( x2 - 1 )$d,x + i(l - y 2 )i>D,y]i>s,x + (1 - v 2 ){x^d,x - yi>D, y )i>s, v } dy, (89) 



where the integral limits are chosen by requiring that the function A to be regular on the axis of symmetry. But 
ips, y = and $s,x = l/(z 2 ~ 1), then 



A SD = 2 



1 



i 2 - y 2 ) 



y{x 2 - l)tpD,x + x(l - y 2 )-4>D,y dy 



(90) 



Remembering for magnetized black hole k = rh = m/ (3, this integral gives 

2(M/0)(m/0) 



A 



SD 



or using (|54bj) we have 



A S d = 



(M/P) 



(m//3)x + z Q y - p 2 + (z + z f 



(91) 



, {[z Q + (m/P)Wp 2 + lz + (rn/P)} 2 

[4 - (m/P) Wp 2 + (z + z ) 2 l. 

- [zo - (m/p)] VP 2 + [z- (m/(3)} 2 - 2(m/f3) ^P 2 + (^ + ^ ) 2 } • 



From (|80|) - (I82|) . the main physical quantities associated with the system are 



4(M//?)z e 

( P 2 + z 2 y/ 2 



F 



(m//3) 



(M//3)p 2 



Pip 



vV + (™//3) 2 (p 2 + ^) 3/2 
(m//3) (M//3)p 2 



4(M//?)z e 

(P 2 + ^o 2 ) 3/2 [VP 2 + (W/3) 2 ' (P 2 + 4) 3/2 
(M/j3)b^_ A p 



2tt 



{ P 2 + z 2 y/ 2 ' 



(92) 

(93) 

(94) 
(95) 



In order to study the behavior of main physical quantities associated with the system we perform a graphical 
analysis of them. In Figs. [T] and [2] we plot as function of p the energy density e, the azimuthal pressure p v , the 
azimuthal electric current density j, the speed v 2 of counterrrotation and the specific angular momentum h 2 for a 
system consisting of a black hole and a Kuzmin-Chazy-Curzon disk of infinite extension in presence of a magnetic 
field with zq = 1, rh = 0.1, M = 0.2, and for values of magnetic field parameter b — 0, 0.5, 1, and 2. We see that 
energy density is a positive quantity in concordance with the weak energy condition, as well as the stress in azimuthal 
direction (pressure). Eq. (|94l) shows that we have pressure for all values of parameters. Since e + p v > 0, the strong 
energy condition is also satisfied. These properties characterize a distribution of matter with the usual gravitational 
attractive property. On the other hand, the dominant energy condition (v < 1) is not satisfied in the central region 
located between the black hole and the photon radius, but the increase of the magnetic field can make this region 
closest the center of the disk. We also observer that these structures present strong instabilities inside the photonic 
orbit. 



VII. BLACK HOLES SURROUNDED BY FINITE DISKS IN AN MAGNETIC FIELD 

Finite thin disks can be obtained introducing oblate spheroidal coordinates (u,v), which are naturally adapted to 
a disk source. They are related to prolate coordinates by x = —iu, y = v, and k = ia, being a other constant, and 
with the Weyl coordinates (p, z) by 



p 2 =a 2 (u 2 + l)(l-v 2 ), 
z = auv, 



(96a) 
(96b) 
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with u > and — 1 < v < 1. In terms of these coordinates the asymptotically flat general solution of Laplace's 
equation P0|) can be written as 



1p = ~ }^ C 2n P 2n (v)q 2n (l 



where c 2n are constants and 



n=0 



q 2n (u) = i 2n+1 Q 2n (iu) 



The function metric A can be obtained by integrating 



A, u = 
A„ = 





w 2 \ 




-v 2 J 


/ u 2 - 






-v 2 ) 



i{u 2 + l)i; 2 u - u(l - v 2 )i> 2 v - 2v(u 2 + l)ip, u i(), v 



v(u 2 + 1)?P 2 U - v(l - v 2 )i> 2 v + 2u(l - v 2 )V>V> ■ 
For magnetized solutions constructed from the Ernst's method the magnetic potential is given by [4£ 

oo 

C 2r , 



A = -ab (u 2 + 1) £ ^y«2>) K2n(«) " 



(97) 

(98) 

(99a) 
(99b) 

(100) 



where q'2 n ( u )( x ) = d<l2n{u) / du. In the absence of magnetic field 6 = these solutions correspond to the Morgan- 
Morgan solutions [HI and represent the exterior gravitational field produced by a finite disklike source localized at 
z = and radius a. In oblate coordinates the disk is localizated in u = and — 1 < v < 1. Following the reference [14| . 
on crossing the disk the coordinate i> changes sign but does not change in absolute value, whereas u is continuous. In 
consequence, as tj) is an even polynomial function of v, it is continuous across of the disk, but its normal derivate v 
is an odd polynomial function of v, and hence discontinuous across of the disk. The same holds for ip and A. This 
implies, as was discussed to the final section II, that these solutions can be interpreted as the exterior gravitational 
field produced by a finite dislike source immersed in a magnetic field. 
For the first two terms n = and n = 1 we have 



A 



-cocot" 1 ^) - i C2 (3v 2 - 1) [(3u 2 + l)cot _1 (u) - 3tt] 



(101a) 
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16a 2 
3 



4p 2 (9u 2 v 2 -u 2 + v 2 -1) 



+-c 2 u(l - v 2 ){27c 2 u 2 v 2 - 3c 2 u 2 + 21c 2 v 2 - 5c 2 + 8c ) cot _1 (w) 
8 

1 



:(c 2 - c ) 2 In 



u 2 + v 2 
1 + u 2 



-777C 2 (1 - v 2 )(12c 2 v 2 + 27c 2 u 2 v 2 - 3c 2 u 2 
lb 



4c 2 



16c ), 



A = abc (v + l) + -abc 2 v (l - v 2 ) [3u (u 2 + l) cot _1 (u) - 3u 2 - 2] 



(101b) 
(101c) 



and ip is given by (I41a|) . In the absence of magnetic field 6 = and Co = c 2 = M/a, where M is the mass of the 
disk, this metric corresponds to the first Morgan-Morgan disk. So the case 6^0 correspond to the first magnetized 
Morgan-Morgan disk with parameter Co = M/(/3a). 

For simplicity, we shall superpose a black hole and the first Morgan-Morgan disk in a magnetic field. According to 
(|99b[) . the interaction term (|b7| between the disk and the black hole in oblate coordinates is given by 



Aszj = 2 



( u 2 + y 2 J \{ V ( u2 + W-D.u + "(1 - V )i)D,v]i)S,u + (1 - V 2 )(w/jd,u ~ V^D,v)^S,vj dv, (102) 



where the integral limits are chosen by requiring that the function A to be regular on the axis of symmetry. But 
tps.v = and 



k 

a{u 2 + 1) ' 



(103) 
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then 



a sd = 2 ± r ^ 2 +i)*d.«+"? dv . m 



a 



Taking k = rh = m/ (3, this integral for the first magnetized Morgan- Morgan disk gives 

Asd = {3c 2 (l - v 2 ) (ucoir^u) l) + (c 2 - c ) In } , 

and the magnetic potential is given by 



A = b(m/(3){y + 1) + abc (v + 1) + habc 2 v (l - v 2 ) [3u (u 2 + l) cot _1 (u) - 3u 2 - 2] . (106) 
Again ip is given by (|41al) . For the combined system the energy density and azimuthal pressure are given by 

r _3nc I f_ (m/(3) \ [e + c 2 (2 - 3p 2 /a 2 )} 




= 4e ,-. tm \^±^Mm. (los) 

1 ' - ^(m//3) 2 + p 2 / aVl-p 2 /* 2 



When cq = ci we have 



P<p 



V 1 a 2 ^ 4 a 2 VW/3) 2 +PV' 
12(M//3) / 37r (M//3)p 2 (m//3) \ 



and the counterrotating speed 



y 2 = (37r/4)(M/(3)p 2 /a 2 + (m//3)/ ^(m//?) 2 + p 2 

- (37r/4)(M//3)p 2 /a 2 - (m//3)/ ^(m/(3) 2 + p 2 ' 

In order to study the behavior of these physical quantities we also perform a graphical analysis of them. In figure 
[3] we show the energy density e, the azimuthal pressure p v , the speed v 2 of counterrrotation, and the specific angular 
momentum h 2 for a system consisting of a black hole and the first magnetized Morgan- Morgan finite disk in presence 
of a magnetic field with radius a = 1, rh = 0.1, cq = c 2 = 0.1, and for values of magnetic field parameter 6 = 
(dashed curves), 0.5, 1, and 2 (dotted curves), as functions of p. We also plot the counterrotating speed and the 
specific angular momentum for a disk of radius a — 100 and the same values of parameters. We see that the energy 
density is a positive quantity in agreement with the weak energy condition, as well as the stress in azimuthal direction 
(pressure). Eq. (|110[) shows that in the case cq = c 2 always we have pressure. The strong energy condition, e+Ptp > 0, 
is also satisfied. On the other hand, the dominant energy condition (v < 1) is not satisfied everywhere on the disk, 
but because the presence of magnetic field decreases the counterrotating speed of particles in all regions of the disk 
we find that after certain value of the magnetic field only at the central region located between the black hole and 
the inner photonic orbit the particles exceed the speed of light whereas in vacuum there are two photonic orbits, 
one inner and other outer. However, imposing that the outer edge of the disk coincides with the outer light radius 
also obtains only an inner region with tachyonic matter [30j, |31( . Since the function F tends to one when p tends a 
zero, the system presents a similar behavior for all values of the parameters. We also observe that these structures 
are stable against radial perturbation except also inside the inner photonic orbit where there are strong instabilities. 
Note that for disks of radius a >• m/(3 (Figs. EId)-[2Ie)) this region is very close to center of the disk. Thus, for a 
reasonably large radius we can consider that these solutions can model approximately the field of a finite disk around 
black hole immersed in a magnetic field of the form (|106[) . 
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VIII. DISCUSSION 



Two relativistic models of thin disk around static black hole in presence of a magnetic field were presented. The 
first model is based in the Kuzmin-Chazy-Curzon inifnite disk and the other in the first Morgan-Morgan finite disk. 
In the first case the source of the magnetic field can be considered the surface electric current density presents on disk 
whereas in second the source is not planar. Thus, these solutions were interpreted as the exterior gravitational field 
produced by a finite dislike source immersed in a magnetic field. 

In both cases, these structures present always a inner region located between the black hole and the photon radius 
where the dominant energy condition (v < 1) is not satisfied. However, the inclusion of magnetic field makes this 
region closest to the black hole. Similarly, was found that these systems are stable against radial perturbation except 
also inside the inner photonic orbit where there are strong instabilities. 

In the case of finite disks, was observed that the increase of magnetic field allows us to have only a photonic orbit 
whereas in vacuum there are two photonic orbit, one inner and other outer. For a reasonably large radius and after 
certain value of the magnetic field was considered that these solutions can model approximately the field of a finite 
disk around black hole immersed in a magnetic field. 

Finally, in order to construct a such system that satisfies all the energy conditions, a model of a black hole surrounded 
by a disk with an inner edge in presence of magnetic field is being considered. 
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FIG. 1. (a) The energy density e and (6) the azimuthal pressure p v for the system black hole and a Kuzmin-Chazy-Curzon 
infinite disk in a magnetic field with zo = 1, rh = 0.1, M = 0.2, and for values of magnetic field parameter 6 = (dashed 
curves), 0.5, 1, and 2 (dotted curves), as functions of p. (c) The azimuthal electric current density j for 6 = (axis p), 0.5, 1, 
and 2 (top curve) and the same value of other parameters. 




FIG. 2. (a) The speed v 2 of counterrrotation and (d) the specific angular momentum h 2 for the system black hole and a 
Kuzmin-Chazy-Curzon infinite disk in a magnetic field for 6 = (top curves), 0.5, 1, and 2 (bottom curves) and the same value 
of other parameters. 
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FIG. 3. (a) The energy density e, (b) the azimuthal pressure p v , (c) the speed v 2 of counterrrotation, and (d) the specific 
angular momentum h 2 for the system black hole and the first Morgan-Morgan finite disk in a magnetic field with radius a = 1, 
m = 0.1, co = C2 = 0.1, and for values of magnetic field parameter b = (dashed curves), 0.5, 1, and 2 (dotted curves), as 
functions of p. (e) The counterrotating speed and (/) the specific angular momentum for disks of radius a — 100 and the same 
value of parameters. 



